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ABSTRACT 


Radiances obtained from the 11-1 window channel by MRIR scanning 
radiometer of NIMBUS II were investigated in order to determine areas 
of time change of Tos which may be related to effective cloud top level 
of the five-day mean cloud state. A statistical model was developed 
that is useful in obtaining a valid specification for effective cloud 
top change, using a stepwise regression program. 

The parameters AT); AZ, and T, appear to provide the necessary 
Spatial and temporal stability to employ the specification scheme over 
a time span of ten days. These three parameters are all highly descrip- 
tive of AT, in terms of basic atmospheric dynamical relationships. It 


is recommended that a test be made of this hypothesis using a 24-hour 


See il |eal10n pnoceaure. 
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I, INTRODUCTION 


The diagnosis of weather features by sequential analysis of satellite 
maps of infrared radiances sensed during the lifetime of the NIMBUS II 
Satellite, presents opportunities for understanding tropical and sub- 
tropical atmospheres generally, and in particular the environment of the 
Peer nical jet stream. These radiances were obtained by the MRIR 
scanning radiometer of NIMBUS II, which scanned in the 11]-u window 
channel (channel 2) and gave an equivalent blackbody temperature (here- 
after called Th). This radiometer also scanned in the 6.7-u water vapor 
channel (channel 1), which corresponds to the effective water-vapor 
temperature (hereafter called 7) in the upper troposphere. For a more 
detailed discussion of MRIR data, the reader is referred to the NIMBUS 
11 User's Guide (Staff, 1966). 

Grid point maps of the two equivalent blackbody temperature fields 
1, and Ty for the period 15 May through 28 July 1966 were provided in 
the form of five-day means through the kind auspices of the Laboratory 
for Atmospheric and Biological Sciences of NASA at Greenbelt, Maryland. 
The use of five-day mean maps presented consecutive fields with no data 
gaps on the quasi-globe defined between 70S-70N latitude (Fig. 1), 
although all of the transient daily features have undoubtedly been 
smoothed out. 

The purpose of this research was to specify statistically the time 
change of Tos denoted hereafter as AT, especially in data-sparce areas 
over the oceans, The actual value of AT, may be related to the change 
in the effective or mean cloud-top level of the five-day mean cloud 


state at grid point (i,j). In turn the field of AT. (i,j) will be 





statistically related to certain derived products which are descriptive 
of basic atmospheric processes. 

Areas of significant AT, values are known to be related to areas of 
variation in the cirriform cloud-top level and/or variation in depth of 
moist layers associated with updrafts or downdrafts in the layer 6-10 km 
[Raschke and Bandeen (1967), Beran, et al. (1968)]. During the course 
of this investigation it became evident that AT, is a most descriptive 
parameter for specifying ATs; and hopefully is amenable to prognostica- 
tion for periods up to 12 hours in advance of the latest radiometric 
swath. If so, the method described here could be altered to provide a 
12-hour prediction for AT. 

Although radiometric data covered the region /0S-70N in Mercator 
grid coordinates, the area of interest for this study was restricted to 


the limits 14.3 to 56.28N for which aren Fiy 


@ wwe t ue ee ‘s¥ 


1) 
a 
- 


were available, At the same time, these limits provided sufficient 
latitude to accomodate each of the subtropical jet streams considered. 
The longitude span of the subtropical jet area was restricted for the 
nurposes of this study to the range 180W to 120W. This area could 
normally be classified as data-sparce. During the data period of this 
study, the area displayed well defined five-day mean subtropical jet 


coverage. 
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TIT. PROCESSING OF DATA 


The radiometric data were distributed over a Mercator grid point 
cnart having i=1,...,/3 grid points in the longitudinal direction, and 
j=l,...,41 rows of latitude covering the meridional span 70S-70N. The 
Mercator base chart had a superimposed square mesh grid, with a mesh 
interval d equivalent to five degrees of longitude, true at the equator 
(Fig. 1). The radiometric data T(i,j) were subjected to an initial 
Filtering wnich eliminated waves of two gridmesh lengths, but left 
unaffected all waves of length L > &d. 

As a possible dynamic predictor, the five-day mean 500 mb Northern 
Hemisphere non-overlapping contour fields for five-day map periods 
synoptic with the radiometric temperature fields were also obtained in 


4 i nm - tales ma mee ~ de ~ 23 wt - =~ .~s Pe t 
puncned cara format covering tne rong tude described, 


Ca 
(D 
VY) 
aud ¢ 
C) 
Gas 
C= 
(f) 
c 


but limited to 18 gridpoint rows poleward of latitude 14.3N (i.e., to 
56.8N in Fig. 1). It should be noted that the contour map fields were 
sufficiently smooth that preliminary filtering was not considered neces- 
sary for these fields. 

Certain operations are defined as representing the time-mean and 
time-differencing operations on temperature data at each grid point of 
the successive maps. These time centered operations are defined below: 


Temperature-mean considered at time k: 
Temperature-time difference at k: 


AT, (isJ) = Teergui ed) eh Thaqti ed) 


1] 





<— 


Keo 


= radiometric five-day mean temperature centered at time k+1/2 


x-direction grid indicator 


= y-direction grid indicator 
= 1,2,...,14, the time-step indicator (half steps are 
centered at times of individual five-day means, Fig. 2). 
OR Sees <o nean > 
i oo . =e See 
V2. kel We, Ke2 kK=14 I4¥%2 ks IS 


Fig. 2. Illustrating time-sequencing of five-day mean periods, and 


the 


time-index scheme used in centered time differencing. 


A similar procedure was followed to form AZ by subtraction of successive 


five-day mean 500-mb contour fields which are also centered at the half 


time-steps indicated in Fig. 2. 


A computer program was generated to carry out these time-centered 


conipulations 


longitudinal 


oe Lt Ve 


over the meridional range 14.3N-50.8N as weii as over the 


extent of this investigation. This resulted in a set of 


fourteen consecutive maps, each of which consisted of five fields of 


time-means and time-differences of the following parameters: AT, ke? 


er lak: 


Ty bk? AZ. These fields were programmed for production 


both in punched card and print-out forns. 


Fifteen consecutive subtropical jets in the eastern Pacific were 


manually determined using the fifteen original five-day mean 500-mb 


contour maps 


for the periods centered at times t=k+1/2, with k=0,...,14. 


The selection of individual jet-latitude coordinates was made by examining 


the test area 180W-120W for the latitude of the tightest contour gradient 


that could be assigned nearest to the subtropical anticyclone on the 


five-day mean 500-mb chart. 
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These five-day mean subtropical jets lay well within the previously 
described latitude range of 14.3N-56.8N in which range most of the large 
contour-change AZ centers could be located. Furthermore the fifteen 
geostrophic mean jets varied only slightly in time and space, so that an 
average of Six grid points lay both to the north and south sides of the 
jet. In order that the statistical test area be consistently more 
extensive than the jet area, the data set was expanded longitudinaily to 
include two gridpoints west of 180!! and three to the east of 120W, Thus 
the sample size for each of the k=1,...,14 map-differenced times was 


constant at 216 in this investigation. 





IIT. STATISTICAL PROCEDURES 


A. EXPERIMENTS FOR DETERMINING PREDICTION TIME-STEPS 

The Biomedical stepwise regression program, BMDO2R [Dixon, 1966], 
was used in the statistical experiments described in the following pages. 
A brief summary of symbols and equations follows. 

When the statistical specification is made for AT.» the predictors 


are understood to be Xa; Xo Xo Xn defined at mean map time k, and each 


grid-point (i,j) as: 


Xa = 1) gary > Ty kety2) = AT k 
Xo = (Zara 2 giyyo) = bey 
(1) 

Xa = To paryo t+ To vayo) = 2To 

Ra = (17 kety2 * Ty e-aye? = 2" ,k 
avy 4 = To eaty2 7 T2,k-1/2! 

The variables defined in (1) appear in the sets of specification 
equations associated with AT, (k) in the form 

AT, = Ay + A,X, + AoX, + AaXa + AgXy (2) 


where A. are the regression coefficients, selected in a stepwise 
manner by the program BMDO2R. 
In order to estimate the variable AT, centered at any one given 
map-time k from the concurrent predictor set in (1), it was necessary 
to determine how large a sample base of dependent data was needed. This 


question was tentatively answered by performing the following three 





Statistical experiments, each of which tested not only the significance 
of the specifications but also the stability of the regressions for at 
least one map-period. 

}. Each of thirteen time-centered periods (i.e,. at time k=1,2,...,13) 
was used to generate a specification equation (2) which could be used 
with the data of the next period (i.e., at times k=2,3,...,14) leading 

to a set of thirteen verifications. This experiment gave the best indi- 
vidual period specifications and also single-period verifications on the 
test used in the experiments described below. The order of entry of the 
predictors varied less from case to case, implying greater stability in 
this system of deriving the specification equations than in tests 2 or 3 
described below. 

é. ihe seven odd-numbered map-time fields were pooled in order to gener- 
Pome inaliescamnnc)| tescneciticationmequation, The coefficients of thts 
equation were used to operate on the variables Xa: Xo» Xo and Xa of 
equation (2) at the even-numbered time periods, so that, in effect, a 
composite single-period forecast scheme was devised. In this case, the 
Specification equation itself did not afford sufficient explained 
variance of the predictand AT, for a valid comparison with the results 

of tests 1] and 3. The transgenerated equation explained even less 
significance when applied to the composited seven-day test data (at times 
k=2,4,...,14). This method indicates the lack of statistically homogene- 
ous samples when drawn from a meteorological period encompassing 13 
five-day mean maps, or 65 days. 

3, The final test undertaken was to use a pooled data base consisting of 


three consecutive periods (e.g., k=1,2,3) in order to generate a three-map 


ils 





composite specification equation for use in predicting AT. for the single 
following period (e.g., k=4). This set of apparent single-step forecasts 
was verified against the known data for time k=4. 

Similarly, the succeeding three-period pooled data base for times 
centered at k=2,3,4 was used for a composite statistical sample to pro- 
vide a verification at map-time k=5, and likewise on to the final period 
k=14, 

Test 3, which used a data base of pooled data from three consecutive 
maps in order to transgenerate an estimated AT, for the fourth in 
succession actually involved on the average a full two-period verification 
test. This test was simply not as effective in describing explained 
variance on the independent sample as test number 1 described above, 
which involved single-map time-steps. Hence all results to be considered 


in Tables 2 and follawina are hased unon the results of test 1 anly. 





As previously noted these comparative tests were primarily based upon 
the relative stability of the verification procedures carried out upon 
independent data. 

The definitions xy; Xo s X. and Xn serve to simplify the description 
of the results in Table 2(a) below. This table refers to the order of 
entry of the four variables in the specification equation associated with 
AT, of equation (2). 

The determination of the F-statistic upon entry, computed at each 
step of the multiple regression program (BMDO2R), forms the basis for the 
order of entry of the variables x. in Table 2. 

The definition of the F-statistic at the Jth step of entry is defined 


as: 


[remaining % var. after step J] 


F [1,n-J-1] = on % exp. var.,step J] - [cum. % exp. var., a) 3) 
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Here n is the number of grid points in the data sample. Since signifi- 
cance is sought at the 95% level of confidence for the Jth variable 
added, Fy 1s compared with a critical value a suggested by Miller (1962) 
as 

> 7 Fay(pagsry(on-dt1), a = .08 

| 

where P is the number of predictors among Kasse aXy employed in the 


regression. Ihe sequence of BG values is easily deduced from tables 


of F-values at the indicated rejection limits for a/(P-J+1). These 


erikical Fi-values are listed below for the Jth selection step. 


WaeBeE 1. Ine critical FS Statistics at the 95% confidence limit 
for entry of the Jth predictor of a total predictor set P=4 
variables. The sample size is n=216 in each case. The last row 
shows the corresponding critical values when only three predictors 
are considered for entry by the stepwise rearession. 





B. TEST FOR AN OPTIMAL NUMBER OF PREDICTORS 
The stepwise regression technique applied to AT. and its predictor-set 


X .,X, on map-day 1 led to the results shown in Table 2(a). Comparison 


yor: 
of the F 


4 
j Upon entry with the critical values in Table 1 at steps J=1,...,4 
indicates that the fourth predictor entering does so with a very low 
F-value and a resulting low probability of rejecting the null hypothesis 
that A, equals 0 in equation (2). However the first three predictors all 


entered with confidence levels in Table 2(a) well in excess of the 95% 


level of belief. 





Part (b) of Table 2 contains the specification results for AT, at 


2 

map-time |] with the variable 27, excluded from the predictor set by 
means of an option of the BMDO2R program. The predictor set in Table 2(b) 
is now reduced to AT,, AZ and aT... An objective rationale for the 
statistical preference of 27, to 2T is given in more detail in IV(C), p 32. 

In Table 2(c), the stability of the specification equation at time k=] 
is tested by correlating the observed AT. values of sample map 2 against 
the values estimated from the transgenerated equation for AT, (k=2), made 


by using the day 1 coefficient matrix listed in Table 2(a). Symbolically 


this operation is performed using the estimator matrix equation: 


van 
AT, (2) = (A,, A,, Aj, As, A 


Oe el hiss Ved oa 4) 


] 
2 


es os 


iJ 62 
Here coefficients A. are determined form the program BMDO2R and are 
listed in column 6 of Table 2(a). 
Some question exists regarding the relative stability of the four- 
predictor equation (4) compared to the three-predictor equation inherently 
given in Table 2(b). This three-predictor transgenerated estimator is 


denoted by the symbol AT, (2) for day 2 and has the form: 


LU (2) =o (CARIES iy 


Lier ae ie 


Here the B. are the coefficients listed alongside the predictors in 


atlies2 (bp), 





(Pebeecwmeopecii 1Cation Of tne predictors anda their order of entry 
in equation (2) and for the map period corresponding to time k=]. 
Part (a) refers to the four predictor specification of AT9; part (b), 
to the three predictor specification of AT». Part (c) refers to the 
verification of Aiz from equation (4) applied to test time k=2. 

Part (d) refers to the verification of Al» from equation (5) applied 
to time k=2. 





















Coeff. of 
Jth var. 


ie. (CORK 
Coeff., 
step J 


OIE 324.625 I) abiiesy) 
0.8149 Bangi 0.02242 


0.8322 19.600 -0.1427) 
J=4 073325 Oe al (4 -0.01269 


~ S ‘ = “sh each 8 ‘ i 
Gee CUM Gren. | ete? Oe at oS Paeuncoe: | 


(b) 324.625 eZ Ay. 60.27 


Se Ji) 0.02214 66.4] 
I Za0z2 -0.04019 ear | 
Constant Term -W---n-c nnn nnn ener nee n ene 


22.1258 
| . | 
(c) 0.6543 160.212 0.93228 | 42.81 
(d) AT, 0.6995 205.054 0.95789 | 48.93 


Variable 
entered 





Fy upon 





Sumy 






21, 
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The verification of a prediction on an independent sample-set such 
as that of map-time 2, is judged to be significant when a large Far 
Statistic results for either iT, On AT, compared to a: the critical 
value necessary at the 95% level (Table 1, with J=1). In other words, 
the question to be answered is whether the verifying F-statistic used 
upon independent test-data is greater than the critical F required for 
the 95% confidence level, i.e., is F(1,216) eZ ) with probability 
}-a > 0.95? If the indicated F-statistic is sufficiently large, then 
there is a probability p > 95% of a linear relationship between AT. and 
the test variables (4) or (5). In parts (c) and (d) of Table 2, both 
F-statistics are sufficiently high to refute the probability of a chance 


relationship, although aT, (based upon equation (5)) is found to have a 


A 
higher F-value than that of AT, of (4), 


C. REQUIREMENTS FOR A USEFUL SPECIFICATION FOR AT, 
A useful specification will be defined as one satisfying the two 
criteria listed below. 
In the first place, the specification equation (4) or (5) at each 


step J must be considered significant at the 95% confidence level using 


the F-statistic upon specification: 


a pees % mean sq. var. expl. by (4) or (5 (6) 


% mean sq. var. unexpl. after use of (4) or (5) 


The F-statistic (J=1,...,P) of (6) is very nearly identical to the 
computation of RS, where Ry is the multiple correlation coefficient at 

the Jth step of entry in Table 2(a,b). The F-statistic of (6) is signifi- 
cant at a confidence level in excess of 95% whenever the F-statistic at 


J=P exceeds the fiducial limit F : 
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F(Pn-P-1) > F (P,n-P-1) , a =.05 (7) 


This requirement is conservatively satisfied at the final step J=P with 
n=216, whenever Rp > 0.50. 

The second criterion for the test variables iT, and/or AT, must 
also satisfy the inequality (7) with P=1, for there now exists only one 
way of applying the specification equation to the verification test data 
of time k+l]. This latter criterion makes a requirement upon the time 
Stability of the specification equation, a concept already used in III(A) 
in deciding upon the sampling procedures to be used. 

The second criterion may be most conveniently judged by computing 
the simple correlation coefficients R(AT, 547) pas or R(AT 51 5) pays 
depending upon whether the four- or three-predictor regression (4) or (5) 
is being tested. The statistics are shown, for example in column 4 of 
Tabie Z (c,d) as they pertain to map 2. 
D. ASSESSMENT OF THE USEFULNESS OF THREE-PREDICTOR RELATIVE 

TO FOUR-PREDICTOR CASES 

The test used in criterion 2 of III(B) was made at every verification- 
time k=2,...,14 with the appropriate coefficient matrices of (4) and (5), 
respectively. It was found that the three-predictor equation led to the 
majority of stable verifications based upon consistently larger values 
of the test-statistics R(AT, »ATo) compared to R(AT, sATo) liisaas 
equivalent to saying that the three-predictor equation (5) gave a 
consistently higher F-level upon verification than did the four-predictor 
equation (4). For the reasons just cited and for economy of display, 
only the tabulations of specification and prediction of aT, (the three- 


predictor mode) will be given in IV, which follows. Furthermore, the 


useful rule for a satisfactory specification, R 2 0.5 will be adhered to, 
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at least in the sense of providing a guide line to a Statistically 


Significant prediction. 
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IV. STATISTICAL RESULTS AND INTERPRETATIONS 


A. TABULAR LISTINGS OF THE INDIVIDUAL SPECIFICATIONS AND VERIFICATIONS 

The results of the specification of AT. for map 1 and the application 
to the prediction of map 2 have already been listed in Table 2. Hence- 
forth the tabular listings will be based upon the use of the three-predictor 
equations for both specification and verification. As an example, the 
results of the specification-prediction process for the case of map 2 


proceeding to map 3 are listed in fable 3. 


TABLE 3. (a) Specification of AT? by means of three predictors and 
their order of entry in equation 2) for the map corresponding to time 
k=2. (b) the predictand verification for AT9 at time k=3, tested using 


trancnaonorvatad antuatinan (8) with the canoffiretante vociultina fram man 2 
We we me ee ee Pe GaSe SS O90 0 t~/ es bh wee wee ww we Ne te lt ee ee ee a Sow Aw dy ¢ 6 wots bow yo onal la 


Step Variable | Mult. Corr. F upon Coeff. of | % exp. 
no. entered | Coeff., ae Jth var. var. 
Step J y after J=3 


0.6480 oa 920 1525586 41.99 
Oe 705 32.748 0.01502 GOR 


J=3 | 0.7140 5.466 -0.03471 3) OFS he) 


Constant Term ----------------------- eee Ne267230 


~~ 
(b) AT. ~ 0.6049 123.476 1.01990 36.59 


wo na nn en en ne en ne nee eee 0.54359 
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TABLE 4. (a) Specification of AT, by means of three predictors and 

their order of entry in equation (2) for the map corresponding to time 

k=3. (b) The predictand verification for ATo at time k=4, tested using the 
transgenerated equation (5) with the coefficients resulting from map 3. 


Step | Variable] Mult. Corr. F upon Coeff. of | % exp. 
| i@e entered | Coeff., ante vibe var. var. 
step J y after J=3 


(a) 4 0.5846 111.1034 I3096/ 34.17 
0.6514 $0. 5562 OF01701 42.43 
J=3 0.6626 By WS ya.e 0.03524 43.90 





Constant Term --------------------+----- -~19.08235 
(b) | J=1 AT, 0.6664 Wgeoee yl ee 1ee58 44.4) 
Constant Term ------------------------- -0.11233 


TABLE 5. (a) Specification of ATp by means of three predictors and their 
order of entry in equation (2) for the map corresponding to time k=4. 

(b) The predictand verification for Ato at time k=5, tested using the 
transgenerated equation (5) with the coefficients resulting from map 4. 















Variable| Mult. Corr. F upon Coeff. of | % exp. 






entered oo entry ee Vai 
150.3109 1.63211 41.26 
027057 34 .8372 0.03241 44.25 
O57 151 -0.04216 
Constant Term maar cea aaa 23.47249 
(oyeed= AT, 109.0053 | 0.83761 | 33.76 


Constant Term ---<-------------<---<--- -0.22636 





TABLE 6. (a) Specification of AT» by means of three predictors and their 
order of entry in equation (2) for the map corresponding to time k=5. 

(b) The predictand verification for AT at time k=6, tested using the 
transgenerated equation (5) with the coefficients resulting from map 5. 


Secpoivotiaple eMule. Cony.)| Fy upen 
no. entered | Coeff., 
step J enury 


0.4798 64.0106 lie2o0o6 Oe) UZ 





(a) 


oo) Fes ds) 0.04833 38.04 
Ob Ie 0.3010 | -0.00859 Sisal hs 


oO e eee oo oe eee ee A lige 















105.8482 0.72407 oe40g 


Constant Term ----------------n en eee ene 0.14901 


TABLE 7. (a) Specification of Al? by means o+ three predictors and their 
order of entry in equation (2) for the map corresponding to time k=6. 

(b) The predictand verification for Al? at time k=7, tested using the 
transgenerated equation (5) with the coefficients resulting from map 6. 

















Step |Variable} Mult. Corr. Fy upon Coeff. of | % exp. 
no. entered Caen. bach Jth Var. var. 
Step J y after J=3 


98.5142 Oe Si25 | Sil oe 
J=2 T, 0.6417 34.9574 } -0.06165 41.18 


0.6666 24533 0.01840 44.44 


Constant Term ------------------------- 33. 55210 
(b) | Je] aT, 0.4619 58.0505 | 0.87025 21.34 


Constant Term ------<---n- een en en ne ne- Oral V2 
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TABLE 8. (a) Specification of ATo by means of three predictors and their 
order of entry in equation (2) for the map corresponding to time k=7. 
(b) The predictand verification for ATo at time k=8, tested using the 
transgenerated equation (5) with the coefficients resulting from map 7. 








Variable} Mult. Corr. | F upon 


Step Coeff. of | % exp. 
no. entered | Coeff., VEnevar var. 
step J SN after J=3 


0.4481 53.7527 0.95690 20.08 
0.5919 49.0300 0.05084 35.03 
Ty 0.6090 6.9456 0.04418 37.09 





J=3 
Constant Term ------------------------ -23.46455 

(b) | J=1 AT 0.6588 164.0839 1.24020 43.40 
Constant Term -----------------.------- -2 .38063 

Weer 2. (a8) Specijicaticn o7 Alm by means of three predictors and their 


order of entry in equation (2) for the map corresponding to time k=8. 
(b) The predictand verification for AT at time k=9, tested using the 
transgenerated equation (5) with the coefficients resulting from map 8. 


Peer. OF | *% Exp. 
Junnevay . Naicn 
after J=3 


1.26707 Sieg? 


0.05761 44.12 





J=3 iy : 0.02007 | 44.44 
Constant Term ------------------------ -12.2448] 
ti) |) set Ne 0.5439 29.59 
| Constant Term ------------------------ 0.75761 
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TABLE 10. (a) Specification of AT, by means of three predictors and their 
order of entry in equation (2) for the map corresponding to time k=9. 

(b) The predictand verification for ATs at time k=10, tested using the 
transgenerated equation (5) with the coefficients resulting from map 9. 












Sele Coen = 201 
i- entry | ster 23 
: 
AT, 0.5368 41.6514 0.69409 28.81 
J=3 T, 0.5785 14.8342 0.06375 33.47 
Constant Term ------------------------ -34 .89899 
(b) 1 | AT. 9.1449 2.10 
Constant Term ------------------------ iesaeo 


eee tee Ye Smeesticatian of Ala shvamoonc of thueemevedictersvand their 
order of entry in equation (2) for the map corresponding to time k=10. 
(b) The predicted verification for AT? at time k=11, tested using the 
transgenerated equation (5) with the coefficients resulting from map 10. 











Varnaple) Mult Corr. Fy upon Coeff. of | % exp. 


Coen i... Jth var. | var. 
step J alo after J=3 


0.3431 28.5494 -0.09331 herd. 
0.3480 0.8302 0.06180 Weld! 


G73306 0.4407 0.00469 | Zou 







Constant Term ------------------------ 52 .30643 


(b) | J=1 aT, 0.1015 


Constant Term ------------------------ 1.17698 


-0.10086 
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TABLE 12. (a) Specification of AT» by means of three predictors and their 
order of entry in equation (2) for the map corresponding to time k=11. 
(b) The predictand verification for ATo at time k=12, tested using the 
transgenerated equation (5) with the coefficients resulting from map 11. 


Step | Variable | Mult. Corr. | F, upon 
no. entered =e anon 
far | at, | o.asce ORCC 0.96817 20.81 
en Ty 0.5057 -0.08226 | 25.57 
oa 


J=| 
J=2 
J=3 OReyaiefe' 7-318 W023/ / Conv 








Coe i. oF 
Jth var. 
after J=3 


% @XPp. 
vars 







(a) 


Constant Term ------------------------ 43 .16838 
(b) | J=1 AT, 0.5698 102.8863 1.14851 32.47 
Constant Term ------------------------ 4.59990 


TABLE 13. (a) Specification of AT2 by means of three predictors and their 
order of entry in equation (2) tor the map Corresponding to Lime kre, 
(b) The predictand verification for AT at time k=13, tested using the 
transgenerated equation (5) with the coefficients resulting from map 12. 


F, upon Coeff. of |} % exp. 
Jth var. var. 
after J=3 


J 
entry 





(a) 0.6457 153.0436 igo OO 41.70 
J=2 AZ 0.7031 Se 2009 0.02136 49.44 


JE Ty 0.7114 5.0147 0.04039 50.61 
Constant Term -------------<---+-------- -20.36539 
(b) a 0.6688 iWemeics 0.76579 | 44.73 





Constant Term -------------n---e---n-- -1.61974 
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TABLE 14. (a) Specification of AT> by means of three predictors and their 
order of entry in equation (2) for the map corresponding to time k=13. 

(b) The predictand verification for AT 9 at time k=14, tested using the 
transgenerated equation (5) with the coefficients resulting from map 13. 


step | Variable | Mult. Corr. F upon CGetieon |. 2exns 
no. entered | Coeff., aes Utne Val’s var. 
step J y after J=3 
(a) AT, 0.5876 112.8828 1.00762 34.53 
re On /ar2 Geel? 0.04857 54.34 
pas | 0.7399 eleZ G202137 54.75 


Constant Term ------------------------ -11.92155 


0.6030 2252797 leeosG 36:30 


Constant Term -----------e-e----------- 1.40262 


oe 












TABLE 15. (a) Specification of AT? by means of three predictors and their 
order of entry in euualion (2) for vac, MAP Bo meeocnomeG) von vile “=e 

(b) The predictand verification for Ailz at time k=], tested using the 
transgenerated equation (5) with the coefficients resulting from map 14. 


Variable | Mult. Corr. Fy upon Coeff. of 
Coeff., Jth var. 


Step J Ee after J=3 


(iS oAo U8, 1.12036 


32.4248 Y203933 


0.6485 22.9468 0.09734 


Constant Term -52.64392 


ad 
AT, Oo] 20628791 Oecroy 


Constant Term -0.5772] 
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In addition it is desirable to list all corresponding sets of AT, 
Specification and verification for each of the twelve subsequent pairs 
of days. For example, Table 4 is analogous to Table 3, except that the 
initial specification is for time k=3 and the verification using the 
transgenerated equation is then applied to k=4. The remaining Tables 
5-15 are listed in sequence in view of a similarity in the pattern of 
the statistics presented in the listings, when useful specification- 
predictions occur. In contrast, a distinct change in this pattern occurs 
when a specification-prediction fails, an example of which occurs at 
map-time k=10. 

In order to exploit the data for maximum usefulness, the specifica- 
tion from k=14 was used as a predictor for k=], the results of which are 


Shown in Table 15. 


B. A SUMMARY OF THE MEAN SPECIFICATILON-PREDICTLON RESULTS 

A statistical parameter that clearly summarizes the specification- 
verification results for the two-map statistical sequences is the final 
percentage explained variance. Hence it is advantageous to list this 
statistic for each specification map together with its corresponding 
Vemimicadt iON. sinis 1S done im Table 16. 

The parameter R2 is identified with the mean overall explained 
variances in the set of twelve useful specifications and corresponding 
verifications, which appear at the end of columns 2 and 7. Then R, the 
mean coefficient of determination, has the value R = 0.6911 and .6082 at 
the indicated stages, respectively. 

Table 16 also lists in column 3, for each specification map k, the 


percentage unexplained variance, which is given by aR 
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TABLE 16. Summary of percentage explained variances of AT 




















at the specification step k, and at step (k+1) in the é 
verification process. 
Map~| 2 Exp. |% Unexp.| Standard PEM alte: 
time! Var. ATo Var. for| Deviation ATo on 
k on map k map k map k+] 
1 68.24) Sileaho 48.93 
2 50.59 | 49.41 36.59 
3 43.90 | 56.10 44.4) 
4 : 51.14 48.86 33.76 
8 | 38.13. 61.87 | 5.82 33.09 
6 Tay 44.44 55.56 21.34 
| 7 + 37,09 G2? 91 | A234 
8 | 44.44 55.56 | 6.62934 | 4.3862 | 9 29,59 
oF /aaes 10 ZG 


9 | 33.47 ° 66.53 
5.14620 
1128.07. =» 71.93 ‘| 5.15380 


2 90.16 49.84 eeGozue 















94.79 Usaha) SalZoze 3.4371 14 oreo) 
14 42.06 57.94 62329772 4.8106 ] Sia Jee 


Mean: 49.11 Oa 89 oe 6710 wai 2171 HOS eS, 
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In addition, the table lists the standard deviation o of AT, 
and its standard error (9) py) for each time k=1,...,14,. The standard 


errors in each case may be computed from equation (8): 


standard deviation 


2 
= : ( standard error \° = 1-R¢ (3) 


It should be noted that the two "unreliable" prediction cases 
(10 and 11), which occur because of a single anomalous specification 


at k=10, have not been included for averaging purposes in Table 16. 


C. NORMAL PATTERNS OF PREDICTOR-ENTRY FOR USEFUL SPECIFICATIONS OF AT, 
The statistical feature most apparent among the useful (R > 0.5) 
Specification cases listed in Tables 2 through 15 was the order of entry 
of the predictors. When the multiple correlation coefficient had values 

in excess of entry of C.e, the order of entry in these tables was 
consistently AT, AZ, and To. As mentioned in IIIA, it is the value 
of the F-statistic upon entry which determines the order of entry, and 
in ten of the twelve cases considered to be useful, the predictors 
entered in this particular order. The other two cases of useful speci- 
fication corresponded to the data samples for maps 6 and 11, in both of 
which the order of entry was AT, ; Tos and AZ. For both of these cases 
the explained variance for specification was acceptable due to the large 
F upon entry at the first step, i.e., that which introduced the variable 
Al, . 
An "unreliable" assessment was given to the specification equation 


for map 9 (Table 10) for which the irregular order of entry (AZ, AT), T5) 


and associated relatively small F upon entry values led to coefficients 


S2 





which could not be meshed with those entering the specification equation 
for map 10, Thus map 10 gave poor verification of map 9, and also poor 
Specification when introduced for the verification of map 11. The 
instability of these two map-sets may be seen at once by reviewing 
Tables 10 and 11, with particular attention given to the coefficient 
matrix of the variables AT, » Wks Ty. Note particularly the systematic 
change in the matrices shown relative to those entered immediately before 
and after maps 9 and 10. 

Using percentage explained variance as a guide to the verification 
of maps 10 and 11, it is apparent from Table 16 that the map-sequence 
10-11 should be deleted as useful specification-prediction combinations. 
When this has been done, the mean fractional explained variance upon 
specification and verification are obtained: 0.4511] and 0.3699, 
ecmeauwwee) Thace méamemarco eben Wie sto lcramce cf vaecentelblofraction- 
al explained variance at both specification and verification, namely 0.25. 

At this point a brief discussion of the option of selecting T, rather 
than Ty and Ty is applicable. The primary reason for this decision was 
the high value of the linear correlation coefficient, r(T, 51); in all 
useful specifications. As an example, the correlation matrix based upon 
map time 2 1s shown in Table 17 below. 


TABLE 17. The predictand-predictor 
linear correlation matrix for map 2. 
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Note in particular that r(T,,T,) = 0.768 so that there exists 
considerable redundancy of similar information when both TY and Ts are 


retained in a four-variable specification. The decision to eliminate 


— ——— 


Ty rather than Ty 


upon the small standard deviation of the observed T > relative to that 


by the control-delete option of BMDO2R is then based 


of T, (approximately in the ratio of 1 to 3). Furthermore, using the 
CONDEL option in testing the three predictor equation, which included 
T, rather than the deleted To resulted in a slight shrinkage of the 


multiple correlation coefficients upon verification. 
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V. EXTENSION OF THE RESULTS TO 24-HOUR SPECIFICATIONS 
A. RE-EXAMINATION OF THE FIVE-DAY SPECIFICATION 
PROCEDURE FOR POSSIBLE EXTENSION 
In the preceding discussion, AT. values centered at the Ti ve-day 
ean map time k were found to be linearly related to the independent 
variables AT, ; AZ and Ty by means of a statistical equation of the 


following form: 


AT. = Bo 1: BAT, i BA AZ 2 BT, (9) 


where each of the variables of the right side is likewise time-centered 
at map time k. In addition, the order of importance of the variables on 
the right side of (9), in thirteen of the fourteen cases, was significant 
using the step-wise regression of the BMDO2R program. The ardar wac 
consistently AT,, AZ, To. 

In order to further establish the Stability of (9), a verification 
test was conducted, based upon the results at time k, to the Five-day 
mean data at time k+l. This resulted in only a slight shrinkage in 
percentage explained variance using the same coefficient matrix 
(By Bs Bas B.) on the independent data of map k+1. 

This stability test was desirable to assess the usefulness of the 
diagnostic equation (9) in deriving a valid Specification. In the 
process of testing (9) for statistical Stability, unreliable specifica- 
tion in one case, k=10, led to two cases of inadequate verification, at 
times k=10 and 11. Statistically this has been explained in IV as due 
to lack of homogeneity in the Five-day mean coefficients By» By Bo Ba, 
which occurred at map-time k=10, relative to the longer-term sample 


means of these coefficients over the full test-period of 14 maps. 
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A pertinent conclusion to be drawn, then, is that the statistical 
model evolved in IV is useful in obtaining a valid specification for 
the effective cloud-top level change for a five-day mean cloud state at 
grid points (i,j), using the Standard Atmosphere lapse rate of aetna 
Thus a negative change in AT, of 6.5K is equivalent to an effective cloud 
top ascent of | kilometer, if fully overcast conditions exist. 

On the other hand, if the mean cloud cover at an oceanic grid point 
remains constant at C, where C < 1.0, and an observed decrease in Ty 
occurs, then the ascent in the effective cloud tops (which covers only 
the fractional area C) is deductible from 

AT 


Ale, 2 Y= 6.5 Kk 


2 
yc 


This formula follows since the window channel radiometer senses the 
ieee cUrs con. CCCs SUPTAcCo™ ommendture tn CToudePree <Peas, enc 
the effective cloud top level over the fractional area C of cloud-cover. 
Finally, the parameters AT, AZ and Ty appear to provide the necessary 
Spatial and temporal stability to employ the specification scheme over 
time spans extending as long as ten days in most of the cases listed 
abies 2,...,15). 
B. 24-HOUR STATISTICAL SPECIFICATION OF AT, AND THE PHYSICAL 
INTERPRETATION OF THE PREDICTOR VARIABLES 
The problem to be considered now is that of the analogous use of 


Strictly 24 hourly data. It would be possible to use the BMDO2R program 


to obtain a diagnostic equation of the form: 
AT, (24) = Cy + C,AT, + CAZ + QT, (10) 


based entirely upon independent variables observed at 24 hourly intervals, 
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and time-centered twelve hours before the latest data report has been 
received. The assumption is that good daily mapped versions of the Ty 
and 1 fields have been received in suitably calibrated digitized Tpp 
form from the designated processor in Washington (NESS), and communicated 
instantly to the satellite data-user in the field. 

The tests applied in this paper do not include 24-hour best-fit 
specification because AT), AT, AZ and Ty were not available in synop- 
tically time-centered form at appropriate times. However, if useful 
specifications were derived on the basis of five-day averaged data of 
the same form, it 1s expected that the specification equations of at 
least the same or better efficiency should result for 24-hour averaged 


data. This concept will be discussed from the point of view that AT. 


depends upon basic atmospheric physical parameters, which are reflected 


a Ceara 
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sequencing of the physical relationships, and therefore of the specifi- 
cations, should be more precise for shorter periods (24 hours) than for 
longer period specifications, such as five-day periods. During this 
longer period, atmospheric diffusion tends to smear the effects of the 
physical relationships which bear on ATs. 

Equation (10) is therefore proposed as a significant regression to 
be applied to 24-hour mean data fields much as such regressions were 
found to be descriptive for five-day means. Further supporting argu- 
ments below are for a 24-hour specification equation of form (10) based 
upon the physical usefulness of the variables on the right side of (10). 

1. The 24-Hour Time Chanae, AT, 

Raschke and Bandeen (1967) have shown that properly calibrated 


values of 1, (when used in conjunction with T,) yield the mean relative 


cy 





humidity in the tropospheric layer extending between six to ten kilo- 
meters. AT, may therefore be interpreted as either a change in the mean 
relative humidity of this column or as a vertical velocity producing the 
Same humidity effect. 

Beran, et al. (1968) showed that 24-hour positive changes in AT, 
are associated with areas of convergence near the 200 mb jet level, with 
resultant downward forced vertical motion giving reduced mean humidities 
in the upper convergent tropospheric layer. The use of five-day means 
in the present investigation, as contrasted with that of Beran, resulted 
in the loss of some discrimination regarding AT, ; Since the five-day mean 
fields of AT, ; near the mean polar front jet had much smaller gradients 


than were found by Beran et al. on a single-day change basis. 


2. The 24-Hour Height Change at 500 mb, AZ 
oe an ee Anne eldcees ehOO mb avodemelayad hore ¢ings 
they have been found to be useful in analysis of cloud-development seen 
in satellite pictures. In that context, Z is usually presented in its 
scale-pattern decomposition: 


(ee tL 


(11) 
AZ = Az 


Here, Since the long-wave field Zon is temporally conservative, it 


follows that AZ ALZons and therefore, 


AZ = Ze, (t=0) - Zon (t=-24). 


SD 


BiStKIDULIONS Of Zon are known to manifest the short wave 500-mb 


vorticity patterns which are superimposed upon the large-scale steering 


patterns, the latter being the Z., fields at 500 mb. In turn, the AZon 


SR 
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patterns are related to photographic features of cloud analysis, e.g., 
the comma-shaped cloud with maxima of positive vorticity advection. 
Thus AZ patterns occurring over the preceding 24 hours represent changes 
in the relative vorticity patterns over the area and are instrumental 
not only in effecting changes in the photographic depiction of cloud 
Masses or bands but also changes in their vertical extent, reflected by 
parameter To. In summary, the parameter AZ, as defined for the past 24 
hours, should be expressive of a cloud-structure change in AT. resulting 
from the vorticity-advection change field within the grid area. 
3. The Large-Scale Steering Parameter, I, 

Finally, in scale decomposition applications, the 500-mb SR 

field is known to represent the steering component of the 500-mb wind 


field insofar as vorticity changes are concerned. Space-average contour 
eee cee mew ialent 1 Ue cae tele were hor available, however 

aht fi oT j 
for any of the tests conducted in this study. Nevertheless, a time- 
smoothed 500-mb map is very nearly equivalent in its long-wave steering 
characteristic to the space-smoothed Zop-maps which may be derived from 
an individual 500-mb map, at the time-centered date. 

The Northern Hemisphere five-day mean contour heights have been 
added to the data file for use in a BMDO2R program for the same set of 
15 data periods as were originally investigated. This was done in a 
manner exactly analogous to that used in forming fields of Tos which is 
a predictand entering the right side of (9). The five-day Northern 
Hemisphere fields of To: for the 1966 period of NIMBUS II, have been 
intensively studied by Professor F. L. Martin, who has found that they 
afford a striking resemblance to the corresponding five-day mean 500-mb 


contour fields, period-for-period north of 15° latitude. 
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This fact suggested the supplementary test of correlating the 
five-day means of Ty and Z over the test area of this investigation 
(216 grid points), based upon a pooled series of five consecutive maps 
of five-day means of both Ty and Z. The ensuing regression analysis gave 


the result: 


Z(gpm) = 1919.0293 + 6.94777. 


withia linear correlation: R(Z, T,) = .8506. 

This result gives a rationale for employing T, in the correlation 
originally tested in the specification equation (9) where, except for AZ, 
only radiometric variables were entered into the regression equation. 

It thus becomes clear that the five-day mean To field has nearly the same 
smoothed properties as the five-day mean Z-field. 

It is proposed that, for the 24-hour specification, the use of 
the time-centered field of Tos when space-smoothed tv grid points on an 
NMC synoptic-scale chart, reflects primarily the SR-advecting component 
of the 500 mb field. Thus in (10), the space-smoothed Ty field should 
by analogy perform the role of the 500-mb SR field in steering the 
individual disturbances. For this reason the space-smoothed T, (from 
the 12-hour earlier swath) has been included as predictor in the 24-hour 
specification equation (10) for AT». 

Finally, it has been concluded, using both the references and 
tests cited here, including related knowledge from satellite-picture 
interpretations, that the three parameters AT); AZ, and Ty which appear 
in (10), are likely to be descriptive of ATs. It is recommended that a 
test be made of this hypothesis using the 24-hour specification procedure. 

From the preceding discussion, each of the parameters just cited 


is representative of a large-scale dynamical factor affecting cloud-structure 
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in the atmosphere. While there is undoubtedly non-linear interaction 

between the physical parameters proposed, the use of a linear regression 

model tends to separate nonlinear effects into linearized ones. Hence 
the following summary of the statistical parameters included in (10) is 
presented here with the suggested physical interpretation which they 
represent: 

(1) AT, is indicative of vertical motion change over the past 24 hours. 

(2) aZ is indicative of the relative vorticity propagation by steering 
and of development during the steering. 

(3) The space-smoothed field 1, is indicative of the space-smoothed 
500-mb field. It is also indicative of a 24-hour mean radiometric 
estimate of the cloud-top level over oceanic grid points. In this 
later connection, Th (iJ) should also be related to average surplus 
(doficit) of radiative eneray in the mean celumn at (4 i) and there- 
fore to the eddy-generation effect on an existing or incipient 


disturbance (Katayama, 1967). 


4 





VI. PROPOSALS FOR FUTURE STUDY; A SUMMARY 


In order for an equation of form (10) to have predictive value for 
a period of 12 hours over a given test area, it will be necessary to have 
a polar-orbiting satellite such as NIMBUS IV, reporting T, and To data 
in equivalent blackbody form at 12 hourly intervals. In this way, the 
24-hour interval centered 12-hours before receipt of the latest report 
would provide the resulting specification equation, and could then be 
applied to a time period centered now at time t=0, concurrent with 
receipt of the latest radiometric maps provided AT, can be computed. 

This scheme would require the t=0 field of Tos smoothed to synoptic 
scale, as well as prognosticated fields of AT, and AZ for the succeeding 
twelve hours. A 24-hour AZ prognosis is presently available from the 
Primm cme cauation moce!. “Prognostic fields o7 Aly are not presentiy 
available; however it is not unreasonable to suggest that with further 
tests and evaluation using satellite-provided waterevapor channel fields 
Such a prognostic scheme could be developed. 

NIMBUS IV included a temperature-humidity infrared radiometer (THIR) 
experiment, the purpose of which was to measure by channels | and 2 
daytime and nighttime surface and cloud-top temperatures and cloud 
mapping; by channel 1, atmospheric water vapor mapping (NIMBUS IV User's 
Guide, 1970). THIR data are archived in separate daytime and nighttime 
swaths. The problem is then one of obtaining prognosticated 12-hour 


fields of T,, so that T)(12) - T)(-12) = AT is known for use in lone 


If the temperature distribution is initially known, the vertical 
distribution as well as total amount of water vapor in a vertical column 


above 500 mb may be inferred (Raschke and Bandeen, 1967). The inference 
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of water-vapor absorber mass below 500 mb has not been very successful 
for radiances in the 6./y channel because the effective-intensity emission 
in channel 1 arises from a level in the upper troposphere (Shenk and 
Salomonson, 1970), Furthermore, characteristic if configurations provide 
information about the upper-tropospheric wind field. Martin and 
Salomonson (1970) have used radiometric information from NIMBUS II MRIR 
to locate the core of the subtropical jet and deduce maximum wind speed. 
Thus since To 1S Indicative of in-cloud temperature and is subject to 
change with changing cloud-structure much as a satellite photograph 
indicates, the most favorable radiometric field to consider for prognosti- 
cation is that of TY). A successful prediction in the T, field is 
reasonable to hypothesize since 1, is descriptive of the mean water- 
vapor mass above 500 mb; that is, at a level generally above the top of 
Romie onc uce, MTN Sc orecnosti cated 1.-fielc, @ fietd at Al 
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would become available for use with the known prognostic field of AZ, and 
the current To feet Gee amin ATs» the most variable radiometric parameter 
of equation (10), would be amenable to computation by the regression 
method proposed here. 

This avenue of approach seems to offer a fruitful use of MRIR radiances 
as predictors of effective cloud-top change. Furthermore, the result 


would afford a three-dimensional aspect to cloud-distribution forecasting. 
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